A MODEL FOR TWIN KILOHERTZ QUASI-PERIODIC 
OSCILLATIONS IN NEUTRON STAR LOW-MASS X-RAY 

BINARIES 

X.-D. Li^ and C.-M. Zhang^ 
lO ' ^Department of Astronomy, Nanjing University, Nanjing 210093, China; lixd@nju.edu.cn 

o ■ 

"^National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, 

China; zhangcm@bao.ac.cn 



O 
(N 



o 
O 

m 
> 



ABSTRACT 

We suggest a plausible interpretation for the twin kiloHertz quasi-periodic 
'sj" i oscillations (kHz QPOs) in neutron star low-mass X-ray binaries. We identify 

QQ i the upper kHz QPO frequencies to be the rotational frequency and the lower 

O I kHz QPOs the standing kink modes of loop oscillations at the inner edge of the 

i/-^ I accretion disk, respectively. Taking into account the interaction between the 

neutron star magnetic field and the disk, this model naturally relates the twin 
QPO frequencies with the star's spin frequencies. We have applied the model 
to four X-ray sources with kHz QPOs detected simultaneously and known spin 
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^ I frequencies. 
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INTRODUCTION 



The kiloHertz quasi-periodic oscillations (kHz QPOs) have been measured in more than 
twenty neutron star low-mass X-ray binaries (NS LMXBs) in their persistent emission with 
the Rossi X-ray Timing Explorer {RXTE), which offered unique insights into the physics of 
strong gravity and dense matter (see van der Klis 2000, 2005 for reviews). In many cases 
the twin kHz QPOs appear simultaneously and the correlations between the pair frequencies 
have been investigated extensively (e.g., Psaltis et al. 1998, 1999b; Belloni et al. 2005). 
Moreover, their frequencies also follow rather tight correlations with other timing features of 
the X-ray emission (Ford & van der Klis 1998; Psaltis et al. 1999a; Belloni et al. 2002). There 
is currently no consensus as to the origin of these QPOs, nor on what physical parameters 
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determine their frequencies, which have been identified with various characteristic frequencies 
in the inner accretion flow (see e.g. Stella & Vietri 1999; Osherovich & Titarchuk 1999; Lamb 
& Miller 2001; Abramowicz et al. 2003). 

With the discovery of the twin kHz QPOs in the accretion-powered miUisecond pulsar 
SAX J1808.4— 3658, it was found that the frequency separation Ai/ is almost half the spin 

frequency (Wijnands et al. 2003). For other sources with detected spin frequencies from 
the burst oscillations, Az/ are shown to be close to either the spin frequencies or half of them 
(van der Klis 2005 and references therein). These findings seem to hint some underlying 
mechanisms relating Ug to the upper and lower kHz QPO frequencies (z/2 and Ui). However, 
the more detailed measurements show that Au is generally inconsistent with a constant 
value of Vs, but varying with z/2 or ui (van der Klis 2000, 2005 and references therein), which 
cast doubts about the validity of the simple beat-frequency interpretation (Miller, Lamb & 
Psaltis 1998; see Lamb & Miller 2001 for a modified version). Osherovich & Titarchuk (1999) 
suggested the lower kHz QPO frequency z/i to be the Keplerian frequency in the disk and 
the higher kHz QPO frequency z/2 the hybrid between I'l and 2i's. Kluzniak et al. (2004) 
showed that the twin kHz QPOs can be explained by a nonlinear resonance in the epicyclic 
motion in the accretion disks, which can lead to the 3:2 ratio for the two main resonances 
(see also Abramowicz & Kluzniak 2001; Abramowicz et al. 2003). 

In this paper we propose an alternative interpretation for the origin of the twin kHz 
QPOs, by considering the interaction between the neutron star magnetic field and the sur- 
rounding accretion disk. We introduce the model in §2, and present its applications to several 
NS LMXBs with the simultaneously detected kHz QPOs and known spin frequencies in §3. 
The possible physical implications and conclusions are given in §4. 

2. MODEL 

Neutron stars in LMXBs generally accretes via an accretion disk. For most part of the 
disk, the plasma rotates in a Keplerian orbit. Close to the NS surface, the stellar magnetic 
fields begin to truncate the disk and control the motion of the plasma, resulting in a non- 
Keplerian boundary layer lying between the magnetosphere corotating with the star and the 
outer Keplerian disk. We assume that the boundary layer is confined by the inner and outer 
radii, it!in and Rq respectively. As conceivable, the plasma corotates with the magnetosphere 
at i?in, and the plasma's motion begins to deviate from Keplerian rotation and take its 
maximum value at Rq (see Fig. 1). 

As for the construction of the model, we identify the upper kHz QPO frequency 1/2 to be 



-3- 



the rotational frequency at Rq, i.e., 

U2 = u{Ro)=^Uk{Ro), (1) 

where ^'k(-Ro) is the Keplerian rotation frequency at Rq, and < ^ ^ 1. The value of ^ 
depends on the rotational frequency distribution inside the boundary layer. Unfortunately, 
there is no analytic solution to the structure of the boundary layer, and the value of C, 
can be evaluated only through numerical calculations. However, based on the qualitative 
characteristics of the rotation rate in the boundary layer, Campbell (1987) has suggested 
the following form for the disk rotation profile close to the magnetosphere 

u(R) = MR) - MRin)(l - a;in) exp[ ^[^/^-~ ^^ ], (2) 

where c^-'in = z^s/i^K(-Rin)- When R ^ Ri^^, ui^R) — > z/k(-R); when R i?in, z^(-R) i^s- 
Equation (2) gives a reasonable description of t^(-R) close to but our interest here focuses 
on the rotational behavior around Rq, at which Eq. (2) fails to be valid for a sufficiently 
wide range of ujin- Instead, we take a slightly modified version of Eq. (2) to account for disk 
rotation, 

u{R) = MR) - MRin){l - ^in) exp[ ^^f/^"^~/^ ]. (3) 

— LOin) 

As an illustration. Fig. 1 shows v{R)/v]^{Rin) against R/Rm for three selected values of 
LOin — 0.2, 0.5 and 0.8 in sohd curves. The dashed curve corresponds to Keplerian rotation 
for comparison. Although the hypothesized form (3) is phenomenological, it shares the main 
features with the numerically calculated results (e.g., Ghosh & Lamb 1979), and the following 
analysis also indicates that its adequate for the disk rotation. 

From Eq. (3) we can determine the location of Rq using the condition du/dR = at 
Rq, and obtain the values of ^ = z/(i?o)/t'K(-Ro) for different values of the so-called fastness 
parameter cu = fs/^iRo)- Figure 2 shows ^ as a function of uj. 

Before discussing the origin of the lower kHz QPOs, we note that in solar physics 
the QPOs of several minute periods in coronal loops have been detected and successfully 
interpreted as the standing kink modes of MHD waves (Aschwanden et al. 1999; Nakariakov 
et al. 1999). Thus, following the similar ideas in treating the QPOs in the coronal loops, we 
suppose that these MHD waves could also be excited at the inner edge of the accretion disks, 
where the poloidal field lines are dragged along the azimuthal direction due to shear motion 
between the star and the disk, and various types of MHD instabilities take place (Ghosh & 
Lamb 1979). The reconnection of the azimuthal magnetic field fines could result in circular 
flux tubes (loops) in the boundary layer, where plasma is confined along a magnetic field line 
with some cross section. It is well known that the coronal loops may be set into oscillations 
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with various modes leading to brightness oscillations (Roberts 2000 for a review). While 
the similar oscillations may also occur in the accretion disks, generated by MHD turbulence, 
here we focus on the fast kink mode of standing MHD waves. These fast waves arise as a 
free mode only for high density loops (Edwin & Roberts 1982), which may be appropriate 
in accretion disks in NS LMXBs. Assume the loop length being the circumference of the 
magnetosphere, 27iRo, the oscillation frequency is 

^=4^' 

where Ck is the kink speed. The wave number k = Ntt /{2ttRo) = N/[2Rq). The integer 
[N — 1) stands for the node number of the vibration along the tube axis, with N = 1 being 
the principal mode. At the inner edge of the accretion disk, the sound speeds are much 
smaller than the azimuthal Alfven speeds ca</), and the plasma inside a loop is much denser 
than its surroundings for the comparable magnetic field strengths, therefore we have the 
kink speed Ck — -\/2ca0 (Roberts 2000), and in turn the frequency of a standing kink mode 



Nca4> 
2v^7ri?o ' 



(5) 



In the following, wc estimate the azimuthal field strength at Rq. We assume that 
the stellar magnetic field lines are initially dipolar and penetrate the accretion disk. The 
differential rotation between the star and the disk generates the azimuthal field component 
Bfj, from the vertical component (Ghosh & Lamb 1979). According to Wang (1995), if 
the growth of is limited by the diffusive decay produced by the turbulent mixing within 
the disk, B^ is given by^ 

MR) - I'^^^^B^iRl (6) 
where the parameter 7 ~ 1 (Aly 1984; Uzdensky, Konigl, & Litwin 2002). 

We assume that the magnitude of Rq is close to that of the Alfven radius where the 
magnetic energy density equals the total kinetic energy density (Davidson & Ostriker 1973), 
i.e., 

= 2r]fyv^\Ro, (7) 

where p is the mass density, and 77 ~ 1. Combining Eqs. (6) and (7) we get the azimuthal 
Alfven speed at i?o, 

CA4>m^iri"'^^^^^^v^{R,). (8) 
^2 



^ Other processes related to the dissipation of the magnetic field give similar expressions of (Wang 
1995). 
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Furthermore, we suggest the lower kHz QPOs to be the principal fast kink mode of the 
standing MHD waves along the field lines at Rq. Inserting Eq. (8) into Eq. (5) we find 
the frequency to be 

-i = 7(f)^/^^^^^-KW. (9) 
Combining Eqs. (1) and (9) we have 

.^,f )(.-.), 

where a = 7(?7/2)^/^ ~ 1 is taken as a free parameter that absorbs the uncertainties in 
determining Rq and B^{Ro). Since the ^{uj) relation is known from Fig. 2, Eq. (10) suggests 
a unique relation between ui and 1^2 for given values of Ug. 

Obviously, the MHD oscillations represent the continuum, i.e., the oscillation frequency 

is a continuous function of a wave number. We ascribe the lower kHz QPOs to be the 
principal fast kink mode with the loop length of L = 2ttRo due to the following facts. 
Firstly, the global modes with = 1 or 2 are easiest to excite in the loops (Roberts, Edwin, 
& Benz 1984). Secondly, for a standing wave the decay time of the MHD oscillations is 
T oc (L/N)^ (Roberts 2000). Apparently, the oscillations with the longer L and smaller N 
last the longer lifetime, and hence are more easily to be detected. Thus, for the fast kink 
mode at Rq, oscillations along have the maximum loop length and smallest wave number 
when L — 27: Ro and N — 1. There may exist the other oscillation modes (e.g. the fast 
sausage mode of standing fast MHD waves, standing slow MHD waves, and propagating 
MHD waves) in the accretion disks, and some of them could be also shown as some kinds 
of QPOs. However, their oscillation frequencies are either too high or too low compared to 
those of the kHz QPOs in NS LMXBs. Additionally it seems that the kink mode is most 
likely to be seen as a standing wave (Roberts 2000) . 



3. COMPARISON WITH OBSERVATIONS 

We have compared the predicted Ui — 1/2 relation with the observations of kHz QPOs in 
four NS LMXBs (4U1608-52, 4U1636-53, 4U1728-34, and 4U1915-05), in which both the 
spin and twin kHz QPO frequencies have been measured (data were provided by T. Bclloni, 
M. Mendez and D. Psaltis). In Fig. 3 the crosses represent the measurements and the solid 
hues stand for the theoretical relations. The spin frequency and the adopted value of a for 
each source are also displayed. In all cases the value of a is of order of unity as expected. The 
theoretical predictions and the measured kHz QPO data match quite well for 4U 1728—34 
and 4U 1915—05, albeit the approximated consistence is obtained for 4U 1608—52 and 4U 
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1636—53, which may be caused by the fact that the real magnetic field structure in these 
sources may be more complicated than in the simplified model considered here. 

4. DISCUSSION AND CONCLUSIONS 

In this paper we have presented a plausible mechanism for the production of twin kHz 
QPOs in NS LMXBs, invoking the magnetic field-accretion disk interaction. Although NSs 
in LMXBs are thought to have weak magnetic fields (~ 10^~^ G), their infiuence on the disk 
rotation has not been paid much attention in the existed works on this subject. Titarchuk, 
Lapdius, & Muslimov (1998) has already suggested that a shock occurs in the transition layer 
where the Kcplcrian disk adjust to sub-Keplerian flow. The disk can undergo various types 
of oscillations under the influence of the gas, radiation, magnetic pressure and gravitational 
force. As to the origin of the oscillations, the boundary layer in our work is similar to but 
closer to the NS compared with the centrifugal barrier region in Titarchuk et al. (1998). 

One of the twin kHz QPO frequencies is usually interpreted as the Keplerian rotation 
frequency at some preferred radii, most likely the inner edge of the disk (e.g., Miller at al. 
1998; Titarchuk & Osherovich 1999). Wc have followed this idea, but suggested that the real 
(non-Keplerian) rotation at Rq leads to the upper kHz QPOs. Unfortunately, since there 
is no general analytic form for the rotational profile within the boundary layer, we have to 
use a phenomenological description for the rotation of disk plasma around the inner edge^. 
Although it seems to be adequate, its validity and accuracy shoTiId be testified more carefully 
in the future by both observational and theoretical investigations. 

Moreover, we interpret the lower kHz QPOs to be the fast kink modes of MHD oscil- 
lations in loops along the field lines at Rq, shared with the similar physical mechanism 
for coronal loop oscillations. Since is generated by B^ through shear motion between 
the NS and the disk, this naturally links the twin kHz QPOs with the stellar spin. But 
it is distinct from the traditional beat-frequency model by the fact that in our work the 
twin OPOs originate from different physical processes. We note that similar hypothesis was 
discussed by Muslimov (1995), who suggested that the QPOs in LMXBs may be caused by 
the excitation of the so-called nonlinear global Alfven modes in the boundary layer plasma. 
It is interesting that these modes were observed in the studies of an ideal MHD spectrum 
of a toroidal plasma (Goedbloed 1975). In addition, the detailed numerical investigation 
of these modes was performed by Appert et al. (1982), and their existence was confirmed 



^Even if we take ^ = 1, that is, Keplerian rotation at Rq, Eq. (10) still fits the observational data fairly 
well. 
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experimentally by Behn et al. (1984) and Evans et al. (1984). 

Recent analysis by Barret, Olive, & Miller (2005) showed that the quahty factor for 
the lower kHz QPOs in 4U 1636—53 increases with frequency up to a factor of ~ 200 when 
ui ~ 850 Hz, then drops at higher frequencies. A ceiling of the lower kHz QPO frequency 
at 920 Hz is also seen. In the frame of the present work, we ascribe these features to the 
evolution of the twisted field lines. Several theoretical studies of the star-disk interaction 
(e.g. Aly 1984, 1985; Uzdensky et al. 2002) have shown that as a dipole field is twisted 
due to differential rotation, the field lines infiate and effectively open up when a critical 
twist angle is attained (Uzdensky et al. 2002). This limits the azimuthal pitch at Rq, 
|-B</,(i?o)/-Bz(-Ro)| = l{^2 — ^s)/^2, to some critical value, say, 7c, demonstrating that a steady 
state of configuration could be established only if the rotation shear 1/2 — I's is small enough. 
These arguments suggest that there may exist a maximum value of 1/2, beyond which the 
field becomes unstable, resulting in decreasing quahty factor of the lower kHz QPOs at 
higher frequencies and a saturation frequency of ui when most of the field lines become open. 

We have presented a qualitative description of the kHz QPO production mechanism 
and a crude quantitative expression of the kHz QPO frequencies, to interpret the observed 
kHz QPO phenomena. As a preliminary exploration, many physical details have not been 
considered, such as in what condition and how much MHD wave energy is produced to 
account for the observed Fourier power spectrum of kHz QPOs. These should be investigated 
more carefully in the future work. 

This work was supported by NSFC under grant number 10025314 and MSTC under 
grant number NKBRSF G19990754. We are grateful for T. Belloni, M. Mendez and D. 
Psaltis for providing the QPO data, and T. P. Li and P. F. Chen for helpful discussions. 
The authors express thanks to an anonymous referee for the critical comments that greatly 
helped improve the manuscript. 
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Fig. 1. — The solid lines show the dimensionless rotational frequencies in the accretion disk 
for uji-a — 0.2, 0.5, and 0.8 (from bottom to top). The dashed line represents the Keplerian 
rotation. 
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Fig. 3. — Predicted relations for twin kHz QPO frequencies and measured frequencies for 
four NS LMXBs 



